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ABSTRACT 

Galactic double white dwarfs were postulated as a source of confusion limited noise for LISA, 
the future space-based gravitational wave observatory. Until very recently, the Galactic population 
consisted of a relatively well studied disk population, a somewhat studied smaller bulge population 
and a mostly unknown, but potentially large halo population. It has been argued that the halo 
population may produce a signal that is much stronger (factor of ~ 5 in spectral amplitude) than the 
disk population. However, this surprising result was not based on an actual calculation of a halo white 
dwarf population but was derived on (i) the assumption that one can extrapolate the halo population 
properties from those of the disk population and (ii) the postulated (unrealistically) high number of 
white dwarfs in the halo. We perform the first calculation of a halo white dwarf population using 
population synthesis models. Our comparison with the signal arising from double white dwarfs in the 
Galactic disk+bulge clearly shows that it is impossible for the double white dwarf halo signal to exceed 
that of the rest of the Galaxy. Using microlensing results to give an upper limit on the content of 
white dwarfs in the halo (~ 30% baryonic mass in white dwarfs), our predicted halo signal is a factor 
of 10 lower than the disk+bulge signal. Even in the implausible case where all of the baryonic halo 
mass is found in white dwarfs, the halo signal does not become comparable to that of the disk+bulge, 
and thus would still have a negligible effect on the detection of other LISA sources. 
Subject headings: binaries: close — Galaxy: halo — gravitational waves — white dwarfs 



1. INTRODUCTION 



The Laser Interferometer Space Antenna (LISA) is a 
proposed joint ESA/NASA mission that will be the first 
spac e-based gravit ational radiation (GR) detector (see 
e.g., Hughes 2006. and references therein ). It has been 



known for some time (IHils et al. 1990) that Galactic 



double white dwarfs will be a prominent source of GR 
for LISA. Thousands of Galactic double white dwarfs 
are expected to be resolved wel l enough to yield their 
masses and orbital par ameters (Nelemans et al. 2001 
2004| |Ruiter et al.|2007[ ) , which will lead to an improved 
understanding of common envelope evolution scenarios 
and the origin of T ype la Sup e rnova and/or subdwarf 
B star progenitors (|Livio||1989| |Webbink 1984 Iben & 
Tutukov|198"4] |Han eta l. 2003). However, a much larger 
number of double white dwarfs (~ 10 7 ) will be detectable 
within the LISA sensitivity range but will be unresolved. 
In fact, close double white dwarfs are so numerous that 
they are expected to dominate the LISA GR signal at low 
frequencies, 6 their signal rising above that of the instru- 
mental noise level and generating confusion- limited noise 
- a confusion 'foreground'. Thus, in order to attempt to 
uncover sources beneath the confusion noise, one must 
construct the expected total GR signal from these bi- 
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naries and remove it from the LISA data stream. To do 
this, it is of considerable importance to determine a priori 
the characteristics which will set the level of signal (e.g., 
masses, orbital periods, location in the Galaxy, etc.) of 
the double white dwarf population. In addition, the level 
of signal arising from Galactic double white dwarfs can 
be useful in constraini ng the structure and extent of the 
Galactic thick disk (Benacquista & Holley-Bockelmann 



20061. 

The contribution of low- and intermediate-mass extra- 
galactic binaries to the LISA signal was shown to be 
rather in significant in comparison w ith the Galactic pop- 
ulations ( Farmer fc Phinney||2003 1 . The GR signal from 
Galactic and extrag alactic black hole MACHO 7 binaries 
was investigated by |Ioka et al. ( 1999[ ) and it was in gen- 
eral found that black hole binaries at cosmological dis- 
tances will have a higher impact on LISA than will the 
halo black hole MACHOs, although the foreground signal 
from Galactic double white dwarfs will still dominate. 

The LISA GR signal arising from the Galactic popula- 
tion of double wh ite dwarfs has been investigated by sev- 



eral groups (e.g.,|Hils et al 



1998; Hil s fc BendeTP OOO 



Postnov fc Prok horovl 



1 990. | 
elemans et a 



2004 



Benacquista et al.l2 004 Edlund et al. 20051 |Timpano e t| 
al.||2006| |Ruiter et al.||2007p . In most previous calcula- 
tions, the Galactic GR signal was calculated for a single 
component disk with the bulge excluded, unt il recently 



where a bulge component has b een considered (Nelemans 



et al. 



2004 Ruiter et al. 20071. We note that whileTEe 



signal arising from the bulge should not be discounted, 
the Galactic disk is the major contributor to the LISA 
GR signal out of the two populations (see § [3] for fur- 
ther explanation). Up until now, a full calculation of the 

7 Massive Astrophysical Compact Halo Object, or sometimes 
MAssive Compact Halo Object. 
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LISA GR signal due to the halo population of binaries Siegel et al. 2002): 



has never been calculated. Hiscock et al. (20001 esti- 



mated that a Galactic halo white dwarf population will 
produce a GR signal significantly exceeding that of a disk 
population, though it was hypothesized by Nelcmans et 
al. (2001 1 that such a strong GR signal from the halo is 
unlikely given the physical characteristics of a > 10 Gyr 
halo stellar population. We demonstrate that the latter 
is indeed true, and argue against the existence of a strong 
halo GR signal with our detailed calculations § 4. 

In this study, we calculate the LISA GR signal aris- 
ing from white dwarf binaries in the Galactic halo. We 
include in our signal calculation any double white dwarf 
within the LISA sensitivity range: f gI = 0.0001 — 0.1 Hz 
(orbital periods between ~ 5.6 hours and 20 s). We note 
here that this work is a followup study to our previous, 
more detailed study of the LISA signal arising from dou- 
ble white dwarfs in the Milky Way disk and bulge. Thus, 
for a thorough description of population synthesis model- 
ing, signal calculations or binary ev olutionary histories, 
we refer the reader to that study ( Ruiter et al.||2007j 
which we will refer to from now on as RBBLW). The 
main objective of this work is to determine whether or 
not the halo population, omitted in most previous stud- 
ies, will provide a significant contribution to the LISA 
GR signal. Since we do not know the contribution of 
white dwarfs, in number or in mass, to the Galactic halo, 
we use two extreme models to bracket our uncertainties. 
In one model we set the halo contribution to zero (e.g., 
no white dwarfs in the halo and the halo white dwarf 
baryonic mass fraction t]b — 0), while in the other we as- 
sume that all of the baryonic halo mass is in white dwarfs 
(?7b = 1). We then calculate models within these brack- 
ets to compare the various halo realizations to that of 
the Galactic disk+bulge population to determine if (and 
at what point) the GR signal from halo double white 
dwarfs becomes significant compared to that of the rest 
of the Galaxy. We discuss our results in context of re- 
cent microlensing experiments in order to constrain our 
intermediate models, and comment on the prospects for 
future LISA observations in § |4j In § [2] we describe our 
calculations, and in §[3] we present our results of the halo 
gravitational wave signal juxtaposed with the Galactic 
(disk and bulg e) gravitational foregr ound calculations of 
RBBLW and jNelemans etal] ( |2004 1 . 



2. MODEL DESCRIPTION 



The Galactic baryonic hal o mass, potential, and shape 
are all poo r ly constrained ([ Zinn 1985; Saha 1985; Mor- 



rison 



1996| |Majewski||1993[ ~ |Morrison et al.||2003| ). In 



general, halo properties are estimated from star counts 
of thousands of intrinsically bright stars, such as red gi- 
ants, blue horizontal branch stars, and RR Lyrae vari- 
ables, with associated photometric or spectroscopic par- 
allaxes. There is strong evidence that at least some of 
the halo was built by destroying s atellite galaxies, result - 
ing in a "stringy" halo structure ( Johnston et al.f l999). 
Some obvious evidence for this may be found in the Mag- 
ellanic Stre am, a tidal arc of stars stretching over 100 ° 
of the sky flPutman et~"aL"1|1998| |Morrison et al.||2000 l. 
Nonetheless, for this first attempt, we employ the canon- 
ical well-mixed spherically symmetric halo with a total 
baryonic mass of 10 9 M@ and a density profile as follows 
( |Ziim|1985||Morrison|1996||Morrison fc Sarajedini |1996 ' 



Phalo OC (1 + r/a ,halo) 



-3.5 



(1) 



where ao,haio is the scale radius of 3.5 kpc. 

Once we have the mass model, we can populate the 
halo with our binaries. We use the StarTrack population 
synth esis code for single and binary evolution (jBclczynski 
et al. 112008 1 to evolve our halo stellar population. We as- 



sume one metallicity and age for the entire halo, although 
we note that the halo has been observed to be comp rised 
of two distinct components varying in metallicity ( |Car-| 
olio et al |20"07t, an d may very well have a triaxial shape 



(e.g., |Helmi ||2004 ). For the halo population we use an 
evolution model which incorporates: i) low metallicity 
Z = 0.0001, ii) a burst of star formatio n at t = Gyr, 



and H i) is evolved through 13 Gyr (e.g., [Schuster et al. 
20061). Our spatial distribution and evolutionary model 



parameters for the disk+bulge population are described 
in detail in § 2 of RBBLW, but we summarize the dif- 
ferences here: i) the disk and bulge stellar populations 
are evolved with near-solar (Z = 0.02) metallicity, ii) the 
disk has a constant star formation history for 10 Gyr, the 
bulge has a constant star formation history for the first 
Gyr with none thereafter, and Hi) both disk and bulge 
populations are 10 Gyr old. The remaining evolutionary 
parameters for disk, bulge and halo populations are the 
same. 

Once the halo stellar population has been evolved, we 
record the physical properties of the close white dwarf 
binaries (orbital periods < 5.6 hours) and calibrate the 
results. We construct a grid of models in which we con- 
strain the total mass of white dwarf stars in the present 
halo relative to the total baryonic halo mass. We choose 
4 different realizations of the halo, keeping the halo mass 
constant, only varying the parameter which sets the frac- 
tion, by mass, of white dwarf stars (single and binary) 
within it (ijb)- Binarity of 50% is assumed. We choose 
mass fractions of 0% (no white dwarfs in the halo) and 
100% as the extreme cases. For intermediate cases, we 
choose 2 mod els of 15% and 3 % (ba se d on microlens- 
i ng re sults of Alcock et al. (|2000 
(120001, and IBrook et al 



Lasserre et al. 
extract all dou- 
ble white dwarfs with GR frequencies within the LISA 
sensitivity range (10~ 4 — 0.1 Hz) to calculate the spec- 
tral amplitudes, and compare the LISA halo spectra to 
that of the Galactic disk+bulge. We calculate the LISA 



timest ream signals using the approach of Rubbo et al. 
( |2004[ ), which are added together yielding the total ob- 
servatory data stream, which is then Fourier transformed 
in order to produce the frequency domain data (the spec- 
tra presented in § [3|. The actual LISA GR spectra is 
obtained using the Benacquista et al. (20041 simulation 
code. We note that all of our white dwarf binaries within 
the LISA sensitivity range have circular orbits. Eccentric 
white dwarf binaries are expected t o arise from dynami - 
cal interactions in globular clusters (Benacquist a |20 0T|, 

Hia 



where the phase space densities arc much higher than in 
the halo or the disk. These eccentric binaries, however, 
could provide a unique opportun ity for learning abou t 
white dwarf structure with LISA (iWillems et al.|2007b . 



8 |Alcock et al. | ( |2000^ suggest that a 20% (by mass) white dwarf 
halo is consistent with their microlensing results, although likely it 
is an overestimate if compared with chemical evolution models. 
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3. RESULTS 

For r/B = 1 in a 10 9 M halo, we obtain 1.5 x 10 9 white 
dwarfs (single and binary); 500 x 10 6 binary white dwarfs 
and out of these 27.5 x 10 6 LISA binary white dwarfs. 
Only 5.5% of double white dwarfs have periods shorter 
than 5.6 hr. Obviously, the above numbers scale down 
linearly wit h r/B- Chemo-dyna mical simulations of the 
Milky Way ( |Brook et al. 2003[ ) have demonstrated that 
a white dwarf dominated halo that is evolved from a 
white dwarf progenito r-dominated initial mass function 
(Chabrier et al. 19961 at early times would lead to an 
overproduction of carbon and nitrogen when compared 
to observed abundances, so the t/b = 1 model is unreal- 
istic. For a realistic upper limit on the white dwarf halo 
contri bution we choose t]b — 0.3 based on |Lasserre et al.| 
( |2000l their figure 2; also Brook et al. (2003) their §5), 
and this yields 8.3 x 10 6 LISA double white dwarfs in 
the halo. For a Galactic disk+bulge with a total stellar 



mass of 6 x 10 10 M Q for all stellar types (Klypin et al. 



2002J, we predict a total of ~ 1.6 x 10 9 white dwarfs 
550 x 10 6 double white dwarfs, out of which only 8% are 
found within the LISA band: 44.5 x 10 6 (see RBBLW § 
3). 

Note that the halo is presumed to include a specific 
fraction of mass in white dwarfs (t/b), therefore the pre- 
dicted numbers of white dwarfs are a direct result of the 
i) adopted halo mass and ii) calculated (with popula- 
tion synthesis) mass and period distributions for halo 
white dwarf binaries with the assumed binary fraction. 
For the rest of the Galaxy, in addition to the calcula- 
tion of double white dwarf properties, we have computed 
the white dwarf mass fraction with the adopted Galac- 
tic field initial mass function (for details see RBBLW) . In 
other words the disk+bulge model results in a true white 
dwarf formation efficiency per unit mass, while in the 
halo model this efficiency is imposed a priori (through the 
straight forward application of observational constraints, 
e.g., MACHOs). 

The disk and halo double white dwarf populations dif- 
fer significantly in numbers and physical properties, due 
to the different environments under which stellar evolu- 
tion proceeds (metallicity, age, star formation history). 
For example, typical average double white dwarf chirp 
masses (M = (M p M s ) 3/5 / (M p + M s ) 1/5 , where M p and 
M s represent the first formed and second formed white 
dwarf masses, respectively) of halo systems are 0.13 M Q 
as compared to 0.19 M for the disk. Also, there are rel- 
atively few short period double white dwarfs in the halo 
(5.5% vs. 8.0% for the disk+bulge) since this population 
is older and a larger number of short period systems have 
merged. In particular, some double white dwarf systems 
form on rather short orbits (e.g., hybrid white dwarfs 9 
with carbon-oxygen white dwarf companions) , and none 
of these systems are found in our 13 Gyr old halo popu- 
lation (see below). 

In Figure [I] we show the number density (per resolv- 
able frequency bin) of LISA white dwarf binaries as a 
function of GR frequency for both the halo and the com- 
bined disk+bulge population of RBBLW. At nearly all 
frequencies, the disk+bulge population outnumbers the 

9 Carbon-oxygen (CO) white dwarfs with a thick helium enve- 
lope. 



halo by nearly a factor of two (~ 45 vs. ~ 28 million 
LISA binaries). However, there is a relative increase in 
the number of halo systems between ~ 0.0002—0.0004 Hz 
(~ 170 — 80 minute orbital periods). This is attributed 
to the fact that there is a relatively higher number of 
RLOF double white dwarfs with hydrogen white dwarf 
donors in the halo population. Double degenerate bi- 
naries with hydrogen white dwarfs take a long time to 
form; on the order of > 10 9 — 10 10 years as opposed to 
~ 10 8 for other, heavier double WD types (e.g., CO-CO, 
hybrid-CO) descended from more massive progenitors. 
Remnant binaries formed from progenitors with more 
massive stars are more common in the (younger) Galac- 
tic disk. Hydrogen white dwarfs are evolved from binary 
progenitors in which a white dwarf (e.g., carbon-oxygen 
or helium) is feeding from a low-mass main sequence star. 
At some stage during the mass transfer, 10 the main se- 
quence donor becomes depleted of enough mass - to a 
mass below that of the hydrogen-burning limit - such 
that it is no longer capable of fusing hydrogen in its core 
and thus becomes degenerate (a hydrogen white dwarf 
is born; specific details about the evolution of these sys- 
tems can be found in § 3.1 of RBBLW). Because all of 
the binaries in the halo are 13 Gyr old, most pre-stable 
RLOF systems have had time to reach contact, and many 
binaries which have evolved from more massive progen- 
itors have since merged. Double white dwarfs involving 
hydrogen donors make up 75 % of our LISA binary white 
dwarfs in the 13 Gyr old halo, spanning a frequency range 
from ~ 0.0002 - 0.0009 Hz. By contrast, white dwarf bi- 
naries with hydrogen white dwarfs only make up 40 % of 
the slightly younger bulge population (9 — 10 Gyr) and 
only 13 % of the younger disk population. Additionally, 
only the disk population contains close white dwarf bi- 
naries with / gr > 4.4 mHz (orbital periods less than 8 
minutes), the majority of which are CO white dwarfs ac- 
creting fromjielnrrriwhite dwarfs (a class of AM CVn bi- 



Warner (|1995|). Regarding massive white dwarfs, 



only 1.5 % of LISA halo double white dwarfs are CO+CO 
systems, where as this fraction is 12 % for the disk. Out 
of the halo CO+CO systems, ^2/3 will merge within a 
Hubble time and have combined masses above 1.4 M , 
making them p otential double degenerate scenario Ty pe 
la supernovae ( flben k Tutukov||1984| |Webbinkj[l984f . 

To obtain a better idea of the effect that location m the 
Milky Way has on the GR amplitude, in Figure[2]we show 
the number density distribution as a function of distance 
from the Sun for the three halo realizations, as well as the 
disk and bulge populations of RBBLW. It is immediately 
obvious that the halo signal is expected to be less than 
that of the rest of the Milky Way, given the extensive 
yet sparse stellar density distribution with respect to the 
bulge and disk. It becomes clear that even though the 
number of halo double white dwarfs are only about a 
factor of 1.6 in number below the disk+bulge combined, 
the GR amplitude is expected to be much weaker given 
the 1/D dependence (Rubbo et al. 20041, and the fact 
that the physical properties of the double white dwarfs 
(chirp masses, separations) in all three populations are 
not drastically different. 

In Figure [3] we show spectral amplitude vs. GR fre- 
quency of three Galactic halo double white dwarf real- 



Mass accretion rates are ~ 10 



M 



yr 
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izations: 15%, 30% and 100% halo models. In addition 
to our halo signals we show the LISA sensitivity curve, 11 
the median signal arising from the Galactic population 
from RBBLW (see their section 3.4.1), and the signal 
of the double white dwarf foreground of Nclcman s et al.| 
( 2004 ) . The shape of the LISA instrumental noise curve 
is a function of the acceleration noise (from LISA's ac- 
celerometers) , positio n noise and the gra vitational wave 
transfer function (see Larson et al.| ( |2000| for a more de- 
tailed description of how the curve is ca lculated). The 
signal presented in |Nelemans et alT ( 2004 1 is a measure of 
the barycentred double white dwari confusion foreg round 
amplitude h (rath er than spectral amplitude hf, see Tim- 
pano et al. ( 2006[ )), and is artificially truncated beyond 
~ 2 mHz where individual binaries become resolved and 
the signal is no longer confusion-limited in that study, 
where as the signal from RBBLW is shown for a range of 
frequencies. We have scaled the amplitude of Nelemans 
et al.| (|2004|) by y/T ohs x a/3/20, accounting for a 1-yr 



observation time for LISA and signal modulation due to 
the motion of LISA, respectively, to arrive at the root 
spectral density (spectral amplitude) hf. We note that 
for even our most extreme (and unphysical) halo realiza- 
tion of itb = 1. the halo signal is well below the signals 



of both |Nelemans et al.| ()2004|) and RBBLW, as well as 
the LISA sensitivity curve. 



4. DISCUSSION 



The amplitude of the GR si gnal from halo white dw arfs 
was previo usly estimated by Hiscock et al. J 2000D. 



In 



that study, |Hiscock et al | (|2o6op used jHils et al.| " fl990 1 

to estimate the number of disk white dwarfs, and arrived 
at the number (both single and binary) A^isk = 6.5 x 10 8 . 
Hiscock et al. ( |2000[ ) assumed that the properties of the 
halo white dwarfs were the same as those of the disk. 
Next they adopted the number of halo white dwarfs to 
be N ha \ — 2 x 10 11 . This was based on the mass of 
MACHOs (2 x 10 11 halo objects with masses ~ 0.1 - 



1 M derived by the MACHO collaboration (Alcock et 



al. 2000 1), and the assumption that all MACHOs are 
white dwarfs. Note that this results in a halo mass of 



10 11 M(7 



; Hiscock et al. (2000) does not present the 
actual numbers of double white dwarfs that are within 
the LISA band for the halo no r for the disk populations. 



Naturally, Hiscock et al. (2000) obtain, due to a very high 
number of halo white dwarfs, a very strong GR signal 
from the halo population. In fact they estimated that 
the level of signal from the halo could be a factor of ~ 5 
stronger than the one arising from the disk. 

Comparison of predictions clearly shows that our num- 
ber for the entire disk+bulge white dwarf population 



1.6 x 10 9 ) is similar to the disk prediction of 



His- 



cock et al. (20001. However, their number for the entire 
halo white dwarf population is ~ 2 orders of magnitude 
highe r than we can presently use based on star count 



data (Zinn 1985 



Sarajedini 11996 



Saha 



Siege 



1985) lMorrison||1996 |Morrison fc 
et al.|2002[): Even if we place the 



entire baryonic halo mass of 10 9 M in white dwarfs we 
obtain only 1.5 x 10 9 white dwarfs in the halo, while for 
a more realistic halo model (t/b ~ 0.15) our number is 
2.25 x 10 8 (see § 3). Such a large discrepancy in the esti- 

11 Online Sensitivity Curve Generator, based on Larson, Hiscock 
& Hellings, http://www.srl.caltech.edu/~shane/sensitivity/ 



mate of the number of white dwarfs in the halo leads to 
a very different final result. In particular, our halo con- 
tribution is never higher than the Galactic disk+bulge 
signal. This marked difference stems from the fact that 
we employ a much smaller (baryon ic) mass of the halo 
(10 9 M ) than |Hiscock et al.| (pOOO) (10 11 M ). 

The simple fact is that with a Milky Way yirial mass of 
appro ximately 10 12 M ( |Klypin et aT|2002| |Li fc White" 



20071, the maximum baryonic mass m the entire Galaxy 
cannot be more than ~ 8 x 10 10 M without violating 
the strong constraint on set by WMAP3 (|Spergel| 



et al. 20071. Given that the Milky Way disk and bulge 



itself is known to have a mass ~6x 10 10 M m ([Klypin et 



aL||2002 1 , it is extremely unlikely that the baryonic halo 
can be as massive as O(10 n ) M . Even if there were 
significant play in the total baryonic mass of the halo, a 
white dwarf population of O(10 n ) M would have had 
to have lost approximately 10 10 M in gas during the 
planetary nebulae phase, and this much gas is likely to 
have been observed. 

We have calculated the LISA gravitational radiation 
signal predicted to arise from double white dwarfs in the 
Galactic halo. In doing so, we have performed the first 
detailed calculation of the halo double white dwarf popu- 
lation and compared its signal to that of the disk+bulge 
population. Thus, for the first time there is a complete 
model Milky Way Galaxy (disk, bulge, halo) calculated 
sclf-consistcntly with the same binary evolution popula- 
tion synthesis code. The evolutionary calculations were 



done with the population synthesis code Star Track (Bel- 



czynski et al.|2008[ ) , and the GR signal calculations were 
obtained with the detailed LISA simulation code of Bc- 



nacquista et al. ( 2004 1 . It was found that the GR signal 
arising from the halo population is significantly smaller 
than that of the rest of the Galaxy, and will not con- 
tribute substantially to the overall Galactic foreground 
signal. Further, if we use recent microlcnsing results in 
order to constrain the mass of the halo white dwarf pop- 
ulation to 30% of the halo baryonic mass, we predict 
that the GR signal arising from the halo is at the level 
of h { « 7.1 X KT 21 Hz 1 / 2 at 1 mHz (see Figure [3). The 
disk+bulge double white dwarf population generates a 
much larger noise level: hf ~ 10~ 19 Hz 1 / 2 at 1 mHz 
both for the StarTrack disk+bulge population (Ruitcr 
et al. 2007) as well as for the previous prediction for the 
combined Galactic disk+bulge populati on obtained with 



a diff erent population synthesis model (Nelemans et al. 
[2004| . 

Therefore, throughout the low- frequency region where 
the disk+bulge signal is confusion- limited, the halo signal 
is a factor of ~ 10 (1.1 in log compared to the foreground 
of Ruiter et al. (2007)) lower than that of the disk+bulge. 
Since we have used an upper limit on the white dwarf 
contribution in the halo, and the actual halo white d warf 
content is probably not high er than ~ 10 — 20% , (Tis 



serand et al. ([2007]); see also |Torres et al.| ( |2002[ ) for an 
estimate and discussion of the number density of halo 
white dwarfs), we predict that the actual halo signal will 
be more than 10 x lower than that of the disk and bulge 
combined. The reduced number of high-frequency halo 
systems compared with the disk population will result 
in a small number of potentially resolvable systems from 
the halo, since no halo systems are found with GR fre- 
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quencies above 4.4 mHz (see Figure [T|. A number of 
Galactic binaries predicted to be resolved with LISA have 
frequencies a bove this value (log(/)~ —2.35; |Nelemans| 
et al. ( 2004 1) . Even for an unrealistic halo model for 
which 77B = 1, the level of the average halo GR signal 
does not significantly approach that of the disk+bulgc, 
and remains below the LISA sensitivity curve. It is clear 
that LISA's ability to detect other sources will not be 
strongly curtailed by halo double white dwarfs, and that 
the GR signal from white dwarf binaries in the rest of 
the Galaxy (primarily the disk) will still constitute the 
prime limiting confusion foreground for LISA. 
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Fig. 1. — Number density (n = (dN/df);8 f is the size of a resolvable frequency bin for a 1 year observation time, 1/T j, s = 30 nHz) for 
the 100 % model halo population of LISA double white dwarfs, and the entire Milky Way disc+bulge population of LISA double white 
dwarfs presented in Ruiter et al. (2007). On average, the disc+bulge population is a factor of ~ 2 greater in number than the 100 % halo 
population. 
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Fig. 2. — Number of LISA white dwarf binaries as a function of distance from the Sun in bin sizes of 100 pc. The Galactic bulge distance 
is concentrated around 8.5 kpc, and there are no more disc systems beyond ~ 50 kpc. The 100% halo realization is nearly an order of 
magnitude below that of the disc distribution for the potentially strongest (closest) GR sources, until > 10 kpc, where the disc begins to 
fall off more steeply. 
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Fig. 3. — LISA gravitational wave spectra (amplitude densities). Shown are the gravitational wave signal for 15%, 30% and 100% (i.e., 
t/b = 0.15, 0.3 and 1, smoothed over 5000 resolvable frequency bins) realizations of the halo population of LISA double white dwarfs 
computed from our population synthesis models; Galactic double white dwarf signal from Ruiter et al. 2007 (red); smoothed Galactic disc 
foreground from Nelemans et al. (2004) truncated beyond ~ 2 mHz where sources start to become resolved (blue dot-dash); and the LISA 
sensitivity curve for a signal to noise ratio of 1 (dashed line). 



